pH-dependent reaction triggering in PmHMGR crystals for time-resolved crystallography.

Serial crystallography and time-resolved data collection can readily be employed to investigate the catalytic mechanism of Pseudomonas mevalonii 3-hydroxy-3-methylglutaryl (HMG)-coenzyme-A (CoA) reductase (PmHMGR) by changing the environmental conditions in the crystal and so manipulating the reaction rate. This enzyme uses a complex mechanism to convert mevalonate to HMG-CoA using the co-substrate CoA and cofactor NAD+. The multi-step reaction mechanism involves an exchange of bound NAD+ and large conformational changes by a 50-residue subdomain. The enzymatic reaction can be run in both forward and reverse directions in solution and is catalytically active in the crystal for multiple reaction steps. Initially, the enzyme was found to be inactive in the crystal starting with bound mevalonate, CoA, and NAD+. To observe the reaction from this direction, we examined the effects of crystallization buffer constituents and pH on enzyme turnover, discovering a strong inhibition in the crystallization buffer and a controllable increase in enzyme turnover as a function of pH. The inhibition is dependent on ionic concentration of the crystallization precipitant ammonium sulfate but independent of its ionic composition. Crystallographic studies show that the observed inhibition only affects the oxidation of mevalonate but not the subsequent reactions of the intermediate mevaldehyde. Calculations of the pKa values for the enzyme active site residues suggest that the effect of pH on turnover is due to the changing protonation state of His381. We have now exploited the changes in ionic inhibition in combination with the pH-dependent increase in turnover as a novel approach for triggering the PmHMGR reaction in crystals and capturing information about its intermediate states along the reaction pathway.

Computational Study of Base-Catalyzed Thiohemiacetal Decomposition in Pseudomonas mevalonii HMG-CoA Reductase.

Thiohemiacetals are key intermediates in the active sites of many enzymes catalyzing a variety of reactions. In the case of Pseudomonas mevalonii 3-hydroxy-3-methylglutaryl coenzyme A reductase (PmHMGR), this intermediate connects the two hydride transfer steps where a thiohemiacetal is the product of the first hydride transfer and its breakdown forms the substrate of the second one, serving as the intermediate during cofactor exchange. Despite the many examples of thiohemiacetals in a variety of enzymatic reactions, there are few studies that detail their reactivity. Here, we present computational studies on the decomposition of the thiohemiacetal intermediate in PmHMGR using both QM-cluster and QM/MM models. This reaction mechanism involves a proton transfer from the substrate hydroxyl to an anionic Glu83 followed by a C-S bond elongation stabilized by a cationic His381. The reaction provides insight into the varying roles of the residues in the active site that favor this multistep mechanism.

Biocompatible, Multi-Mode, Fluorescent, T2 MRI Contrast Magnetoelectric-Silica Nanoparticles (MagSiNs), for On-Demand Doxorubicin Delivery to Metastatic Cancer Cells.

There is a need to improve current cancer treatment regimens to reduce systemic toxicity, to positively impact the quality-of-life post-treatment. We hypothesized the negation of off-target toxicity of anthracyclines (e.g., Doxorubicin) by delivering Doxorubicin on magneto-electric silica nanoparticles (Dox-MagSiNs) to cancer cells. Dox-MagSiNs were completely biocompatible with all cell types and are therapeutically inert till the release of Doxorubicin from the MagSiNs at the cancer cells location. The MagSiNs themselves are comprised of biocompatible components with a magnetostrictive cobalt ferrite core (4−6 nm) surrounded by a piezoelectric fused silica shell of 1.5 nm to 2 nm thickness. The MagSiNs possess T2-MRI contrast properties on par with RESOVIST™ due to their cobalt ferrite core. Additionally, the silica shell surrounding the core was volume loaded with green or red fluorophores to fluorescently track the MagSiNs in vitro. This makes the MagSiNs a suitable candidate for trackable, drug nanocarriers. We used metastatic triple-negative breast cancer cells (MDAMB231), ovarian cancer cells (A2780), and prostate cancer cells (PC3) as our model cancer cell lines. Human umbilical vein endothelial cells (HUVEC) were used as control cell lines to represent blood-vessel cells that suffer from the systemic toxicity of Doxorubicin. In the presence of an external magnetic field that is 300× times lower than an MRI field, we successfully nanoporated the cancer cells, then triggered the release of 500 nM of doxorubicin from Dox-MagSiNs to successfully kill >50% PC3, >50% A2780 cells, and killed 125% more MDAMB231 cells than free Dox.HCl. In control HUVECs, the Dox-MagSiNs did not nanoporate into the HUVECS and did not exhibited any cytotoxicity at all when there was no triggered release of Dox.HCl. Currently, the major advantages of our approach are, (i) the MagSiNs are biocompatible in vitro and in vivo; (ii) the label-free nanoporation of Dox-MagSiNs into cancer cells and not the model blood vessel cell line; (iii) the complete cancellation of the cytotoxicity of Doxorubicin in the Dox-MagSiNs form; (iv) the clinical impact of such a nanocarrier will be that it will be possible to increase the current upper limit for cumulative-dosages of anthracyclines through multiple dosing, which in turn will improve the anti-cancer efficacy of anthracyclines.

A Quantum-Guided Molecular Mechanics Force Field for the Ferrocene Scaffold.

Ferrocene derivatives have a wide range of applications, including as ligands in asymmetric catalysis, due to their chemical stability, rigid backbone, steric bulk, and ability to encode stereochemical information via planar chirality. Unfortunately, few of the available molecular mechanics force fields incorporate parameters for the accurate study of this important building block. Here, we present a MM3* force field for ferrocenyl ligands, which was generated using the quantum-guided molecular mechanics (Q2MM) method. Detailed validation by comparison to DFT calculations and crystal structures demonstrates the accuracy of the parameters and uncovers the physical origin of deviations through excess energy analysis. Combining the ferrocene force field with a force field for Pd-allyl complexes and comparing the crystal structures shows the compatibility with previously developed MM3* force fields. Finally, the ferrocene force field was combined with a previously published transition-state force field to predict the stereochemical outcomes of the aminations of Pd-allyl complexes with different amines and different chiral ferrocenyl ligands, with an R2 of ∼0.91 over 10 examples.

Automated fitting of transition state force fields for biomolecular simulations.

The generation of surrogate potential energy functions (PEF) that are orders of magnitude faster to compute but as accurate as the underlying training data from high-level electronic structure methods is one of the most promising applications of fitting procedures in chemistry. In previous work, we have shown that transition state force fields (TSFFs), fitted to the functional form of MM3* force fields using the quantum guided molecular mechanics (Q2MM) method, provide an accurate description of transition states that can be used for stereoselectivity predictions of small molecule reactions. Here, we demonstrate the applicability of the method for fit TSFFs to the well-established Amber force field, which could be used for molecular dynamics studies of enzyme reaction. As a case study, the fitting of a TSFF to the second hydride transfer in Pseudomonas mevalonii 3-hydroxy-3-methylglutaryl coenzyme A reductase (PmHMGR) is used. The differences and similarities to fitting of small molecule TSFFs are discussed.

Proofreading experimentally assigned stereochemistry through Q2MM predictions in Pd-catalyzed allylic aminations.

The palladium-catalyzed enantioselective allylic substitution by carbon or nitrogen nucleophiles is a key transformation that is particularly useful for the synthesis of bioactive compounds. Unfortunately, the selection of a suitable ligand/substrate combination often requires significant screening effort. Here, we show that a transition state force field (TSFF) derived by the quantum-guided molecular mechanics (Q2MM) method can be used to rapidly screen ligand/substrate combinations. Testing of this method on 77 literature reactions revealed several cases where the computationally predicted major enantiomer differed from the one reported. Interestingly, experimental follow-up led to a reassignment of the experimentally observed configuration. This result demonstrates the power of mechanistically based methods to predict and, where necessary, correct the stereochemical outcome.

HSP90 inhibitors reduce cholesterol storage in Niemann-Pick type C1 mutant fibroblasts.

Niemann-Pick type C1 (NPC1) disease is a lysosomal lipid storage disorder caused by mutations of the NPC1 gene. More than 300 disease-associated mutations are reported in patients, resulting in abnormal accumulation of unesterified cholesterol, glycosphingolipids, and other lipids in late endosomes and lysosomes (LE/Ly) of many cell types. Previously, we showed that treatment of many different NPC1 mutant fibroblasts with histone deacetylase inhibitors resulted in reduction of cholesterol storage, and we found that this was associated with enhanced exit of the NPC1 protein from the endoplasmic reticulum and delivery to LE/Ly. This suggested that histone deacetylase inhibitors may work through changes in protein chaperones to enhance the folding of NPC1 mutants, allowing them to be delivered to LE/Ly. In this study, we evaluated the effect of several HSP90 inhibitors on NPC1I1061T skin fibroblasts. We found that HSP90 inhibition resulted in clearance of cholesterol from LE/Ly, and this was associated with enhanced delivery of the mutant NPC1I1061T protein to LE/Ly. We also observed that inhibition of HSP90 increased the expression of HSP70, and overexpression of HSP70 also reduced cholesterol storage in NPC1I1061T fibroblasts. However, we did not see correction of cholesterol storage by arimoclomol, a drug that is reported to increase HSP70 expression, at doses up to 0.5 mM. The increase in other chaperones as a consequence of HSP90 improves folding of NPC1 protein and relieves cholesterol accumulation in NPC1 mutant fibroblasts.

Inhibition of Histone Deacetylases 1, 2, and 3 Enhances Clearance of Cholesterol Accumulation in Niemann-Pick C1 Fibroblasts.

Niemann-Pick disease type C1 (NPC1) is a rare genetic cholesterol storage disorder caused by mutations in the NPC1 gene. Mutations in this transmembrane late endosome protein lead to loss of normal cholesterol efflux from late endosomes and lysosomes. It has been shown that broad spectrum histone deacetylase inhibitors (HDACi's) such as Vorinostat correct the cholesterol accumulation phenotype in the majority of NPC1 mutants tested in cultured cells. In order to determine the optimal specificity for HDACi correction of the mutant NPC1s, we screened 76 HDACi's of varying specificity. We tested the ability of these HDACi's to correct the excess accumulation of cholesterol in patient fibroblast cells that homozygously express NPC1 I1061T , the most common mutation. We determined that inhibition of HDACs 1, 2, and 3 is important for correcting the defect, and combined inhibition of all three is needed to achieve the greatest effect, suggesting a need for multiple effects of the HDACi treatments. Identifying the specific HDACs involved in the process of regulating cholesterol trafficking in NPC1 will help to focus the search for more specific druggable targets.

Microsecond timescale MD simulations at the transition state of PmHMGR predict remote allosteric residues.

Understanding the mechanisms of enzymatic catalysis requires a detailed understanding of the complex interplay of structure and dynamics of large systems that is a challenge for both experimental and computational approaches. More importantly, the computational demands of QM/MM simulations mean that the dynamics of the reaction can only be considered on a timescale of nanoseconds even though the conformational changes needed to reach the catalytically active state happen on a much slower timescale. Here we demonstrate an alternative approach that uses transition state force fields (TSFFs) derived by the quantum-guided molecular mechanics (Q2MM) method that provides a consistent treatment of the entire system at the classical molecular mechanics level and allows simulations at the microsecond timescale. Application of this approach to the second hydride transfer transition state of HMG-CoA reductase from Pseudomonas mevalonii (PmHMGR) identified three remote residues, R396, E399 and L407, (15-27 Å away from the active site) that have a remote dynamic effect on enzyme activity. The predictions were subsequently validated experimentally via site-directed mutagenesis. These results show that microsecond timescale MD simulations of transition states are possible and can predict rather than just rationalize remote allosteric residues.

Stereoselectivity Predictions for the Pd-Catalyzed 1,4-Conjugate Addition Using Quantum-Guided Molecular Mechanics.

The conjugate addition of aryl boronic acids to enones is a powerful synthetic tool to introduce quaternary chiral centers, but the experimentally observed stereoselectivities vary widely, and the identification of suitable substrate-ligand combinations requires significant effort. We describe the development and application of a transition-state force field (TSFF) by the quantum-guided molecular mechanics (Q2MM) method that is validated using an automated screen of 9 ligands, 38 aryl boronic acids, and 22 enones, leading to a MUE of 1.8 kJ/mol and a R2 value of 0.877 over 82 examples. A detailed error analysis identified the structural origin for the deviations in the small group of outliers. The TSFF was then used to predict the stereoselectivity for 27 ligands and 59 enones. The vast majority of the virtual screening results are in line with the expected results. Selected results for 6-substituted pyrox ligands, which were not part of the training set, were followed up by density functional theory and experimental studies.

Design, Synthesis, and Evaluation of a Luminescent Cholesterol Mimic.

The development of new chemical tools with improved properties is essential to chemical and cell biology. Of particular interest is the development of mimics of small molecules with important cellular function that allow the direct observation of their trafficking in a cell. To this end, a novel 15-azasterol has been designed and synthesized as a luminescent cholesterol mimic for the monitoring of cholesterol trafficking. The brightness of this probe, which is ∼32-times greater than the widely used dehydroergosterol probe, is combined with resistance to photobleaching in solution and in human fibroblasts and an exceptionally large Stokes-like shift of ∼150-200 nm. The photophysical properties of the probe have been studied experimentally and computationally, suggesting an intersystem crossing to the triplet excited state with subsequent phosphorescent decay. Molecular dynamics simulations show a similar binding mode of cholesterol and the azasterol probe to NPC proteins, demonstrating the structural similarity of the probe to cholesterol.

Transition State Force Field for the Asymmetric Redox-Relay Heck Reaction.

A transition state force field (TSFF) was developed using the quantum-guided molecular mechanics (Q2MM) method to describe the stereodetermining migratory insertion step of the enantioselective redox-relay Heck reaction for a range of multisubstituted alkenes. We show that the TSFF is highly predictive through an external validation of the TSFF against 151 experimentally determined stereoselectivities resulting in an R2 of 0.89 and MUE of 1.8 kJ/mol. In addition, limitations in the underlying force field were identified by comparison of the TSFF results to DFT level calculations. A novel application of the TSFF was demonstrated for 31 cases where the enantiomer predicted by the TSFF differed from the originally published values. Experimental determination of the absolute configuration demonstrated that the computational predictions were accurate, suggesting that TSFFs can be used for the rapid prediction of the absolute stereochemistry for a class of reactions. Finally, a virtual ligand screen was conducted utilizing both the TSFF and a simple molecular correlation method. Both methods were similarly predictive, but the TSFF was able to show greater utility through transferability, speed, and interpretability.

2-Hydroxypropyl-ß-cyclodextrin is the active component in a triple combination formulation for treatment of Niemann-Pick C1 disease.

Niemann-Pick type C1 (NPC1) disease is a fatal neurovisceral disease for which there are no FDA approved treatments, though cyclodextrin (HPßCD) slows disease progression in preclinical models and in an early phase clinical trial. Our goal was to evaluate the mechanism of action of a previously described combination-therapy, Triple Combination Formulation (TCF) - comprised of the histone deacetylase inhibitor (HDACi) vorinostat/HPßCD/PEG - shown to prolong survival in Npc1 mice. In these studies, TCF's benefit was attributed to enhanced vorinostat pharmacokinetics (PK). Here, we show that TCF reduced lipid storage, extended lifespan, and preserved neurological function in Npc1 mice. Unexpectedly, substitution of an inactive analog for vorinostat in TCF revealed similar efficacy. We demonstrate that the efficacy of TCF was attributable to enhanced HPßCD PK and independent of NPC1 protein expression. We conclude that although HDACi effectively reduce cholesterol storage in NPC1-deficient cells, HDACi are ineffective in vivo in Npc1 mice.

Application of Q2MM to predictions in stereoselective synthesis.

Quantum-Guided Molecular Mechanics (Q2MM) can be used to derive transition state force fields (TSFFs) that allow the fast and accurate predictions of stereoselectivity for a wide range of catalytic enantioselective reactions. The basic ideas behind the derivation of TSFFs using Q2MM are discussed and the steps involved in obtaining a TSFF using the Q2MM code, publically available at github.com/q2mm, are shown. The applicability for a range of reactions, including several non-standard applications of Q2MM, is demonstrated. Future developments of the method are also discussed.

pHP-Tethered N-Acyl Carbamate: A Photocage for Nicotinamide.

The synthesis of a new photocaged nicotinamide having an N-acyl carbamate linker and a p-hydroxyphenacyl (pHP) chromophore is described. The photophysical and photochemical studies showed an absorption maximum at λ = 330 nm and a quantum yield for release of 11% that are dependent upon both pH and solvent. While the acyl carbamate releases nicotinamide efficiently, a simpler amide linker was inert to photocleavage. This photocaged nicotinamide has significant advantages with respect to quantum yield, absorbance wavelength, rate of release, and solubility that make it the first practical example of a photocaged amide.

Single-Flask Multicomponent Synthesis of Highly Substituted α-Pyrones via a Sequential Enolate Arylation and Alkenylation Strategy.

Trisubstituted α-pyrones are obtained by a Pd-catalyzed three-component, single-flask operation via an α-arylation, subsequent α-alkenylation, alkene isomerization, and dienolate lactonization. A variety of coupling components under mild conditions afforded isolated yields of up to 93% of the pyrones with complete control of regioselectivity. Metal dependence was noted for three of the steps of the pathway. Utility of the pyrone products was demonstrated by further transformations providing convenient access to polyaromatic compounds, exhibiting broad molecular diversity.

Structural insights into HDAC6 tubulin deacetylation and its selective inhibition.

We report crystal structures of zebrafish histone deacetylase 6 (HDAC6) catalytic domains in tandem or as single domains in complex with the (R) and (S) enantiomers of trichostatin A (TSA) or with the HDAC6-specific inhibitor nexturastat A. The tandem domains formed, together with the inter-domain linker, an ellipsoid-shaped complex with pseudo-twofold symmetry. We identified important active site differences between both catalytic domains and revealed the binding mode of HDAC6 selective inhibitors. HDAC inhibition assays with (R)- and (S)-TSA showed that (R)-TSA was a broad-range inhibitor, whereas (S)-TSA had moderate selectivity for HDAC6. We identified a uniquely positioned α-helix and a flexible tryptophan residue in the loop joining α-helices H20 to H21 as critical for deacetylation of the physiologic substrate tubulin. Using single-molecule measurements and biochemical assays we demonstrated that HDAC6 catalytic domain 2 deacetylated α-tubulin lysine 40 in the lumen of microtubules, but that its preferred substrate was unpolymerized tubulin.

Stereoselectivity in (Acyloxy)borane-Catalyzed Mukaiyama Aldol Reactions.

The origin of diastereo- and enantioselectivity in a Lewis acid-catalyzed Mukaiyama aldol reaction is investigated using a combination of dispersion corrected DFT calculations and transition state force fields (TSFF) developed using the quantum guided molecular mechanics (Q2MM) method. The reaction proceeds via a closed transition structure involving a nontraditional hydrogen bond that is 3.3 kJ/mol lower in energy than the corresponding open transition structure. The correct prediction of the diastereoselectivity of a Mukaiyama aldol reaction catalyzed by the conformationally flexible Yamamoto chiral (acyloxy) borane (CAB) requires extensive conformational sampling at the transition structure, which is achieved using a Q2MM-derived TSFF, followed by DFT calculations of the low energy conformational clusters. Finally, a conceptual model for the rationalization of the observed diastereo- and enantioselectivity of the reaction using a closed transition state model is proposed.